BVRJc photometry of the GRB 980703 and GRB 990123 host 
galaxies 
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Abstract. We present a photometry of GRB 980703 and GRB 990123 host galaxies which was per- 
formed about 20 days and a half year after gamma-ray bursts occured, respectively. The contributions 
of the optical transiets (OT) were negligible in both cases. We derived broad band BVRcIc spectra 
of the host galaxies and compared them to continuum spectra of different Hubble-type galaxies and 
averaged spectra of starburst galaxies. For Hq = 60 km s'^Mpc"^ and three Friedmann cosmolog- 
ical models with matter density and cosmological constant parameters {Qrn,^A) = (1,0),(0,0), (0,1) 
we estimated Ms^^^t and star-forming rates (SFRs) using the fluxes in photometric bands for the host 
galaxies. Within the range of cosmological parameters our estimates of the absolute magnitudes are: 
Mb,,,, = -20.60 ^ -21.73 for the GRB 980703 host galaxy and Mb,,,, = -20.20 ^ -21.82 for the 
GRB 990123 host galaxy. We obtained estimates of K-conection values and absolute magnitudes of the 
host galaxies using SEDs for star-forming galaxies. 

Key words: gamma-rays: bursts: cosmology: observations - photometry: individual (GRB 980703, 
GRB 990123): galaxy - starbursts 



1. Introduction 

The origin of cosmic 7-ray bursts (GRBs) is still one of 
the outstanding problems in modem astronomy. Recently 
the study of GRBs has been revolutionized by the dis- 
covery of X-ray (Costa et al. 1997), optical (van Paradijs 
et al. 1997) and radio transients (Frail et al. 1997). The 
follow-up of optical transients (OT) resulted in the major- 
ity not only with OTs' light curves but also their redshift 
measurements and the detection of GRB hosts. By now 
about 15 optical afterglows have been detected, almost 
all within a fraction of an arc second of very faint galax- 
ies, with typical R-band magnitudes of 22-26 and over. 
About 10 redshifts have been measured in the range from 
z = 0.6 to z = 3.4. The long duration GRBs appear to 
be associated with star forming regions in the galaxies. 
The studies of the host galaxies and their properties ad- 
vance understanding of the nature of the progenitor sys- 
tems. In some proposed models an association of bursts 
with explosions of massive stars become popular (see ref. 
Paczynski, 1999 and MacFadyen A. & Woosley S. E., 
1999 and references therein). 

Here we report our direct BVRcIc imagings of the 



GRB 980703 optical counterpart which were obtained 
after observations of Bloom et al. (1998). The latest ob- 
servations of this object were obtained in the infrared 
JHK bands by Bloom et al. (1998) In the case of the 
GRB 990123 host galaxy our BVRcIc observations are 
the latest. In this paper we compare the BVRcIc spectra 
of the GRB 980703 and GRB 990123 host galaxies with 
spectral energy distributions of normal galaxies of differ- 
ent Hubble types and extend the comparision to averaged 
spectra (SEDs) of S 1 and S2 star-forming galaxies (Con- 
noly et al., 1995). We obtained estimates of K-correction 
values and absolute magnitudes of the host galaxies us- 
ing spectra of normal Hubble-type galaxies and SEDs 
for star-forming galaxies. Rough estimates of low lim- 
its for star-forming rates (SFR) in the GRB 980703 and 
GRB 990123 host galaxies were performed using the 
continuum luminosity at A = 2800A by extrapolating 
for the GRB 990123 host galaxy between Rc and Ic 
bands and by value of the flux in the V band for the 
GRB 980703 host galaxy. 
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Table 1: Photometry of the host galaxies of the GRB 980703 and GRB 990123 



Host 


Date UT 


Band 


Exp. 


Magnitude, obs 


Seeing 








(s) 






GRB 980703 


24.05 Jul. 1998 


B 


480 


23.40 ± 0.12 


1".3 




24.06 Jul. 


V 


320 


22.85 ±0.10 


1".2 




24.06 Jul. 


Rc 


300 


22.44 ± 0.08 


1".2 




24.07 Jul. 


Ic 


360 


22.26 ±0.18 


1".2 


GRB 990123 


8.85 Jul. 


B 


600 


24.90 ±0.16 


1".5 




8.86 Jul. 


V 


600 


24.47 ±0.13 


1".3 




8.84 Jul. 1999 


Rc 


600 


24.47 ±0.14 


I'M 




8.87 Jul. 


Ic 


600 


24.06 ±0.3 


1".3 



2. Observations and data reduction 

Observations of the host galaxies of GRB 980703 
and GRB 990123 were performed using the pri- 
mary focus CCD photometer of the 6m tele- 
scope of SAO RAS. It was carried out with 
the standard (Johnson-Kron-Cousins) photomet- 
ric BVRcIc system. The direct BVRcIc images 
for the GRB 980703 host galaxy are given in 
http : //www. sao . ru/'sokolov/GRB/gSOTOS .html 
Table [l] presents the summary of observations. [[ We 
performed photometric calibrations using the Landolt 
standard fields (Landolt, 1992): PG1633, PG1657, 
PG2331, PG2336 for the host galaxy of GRB980703, 
and PG2213, SAl 10 for the host galaxy of GRB990123. 
To estimate the Milky Way redenning for the host 
galaxy of GRB980703, we used extinction values from 
CardeUi, Clayton & Mathis (1989). We derived 0'".251, 
0".188, 0"'.141 and O^^.OQO in our photometric bands 
BVRcIc respectively. Our direct BVRcIc imagings of 
the GRB 980703 optical counterpart were obtained after 
observations of Bloom et al. (1998) on July 18 UT. The 
most recent observations of this object were obtained 
in the infrared JHK bands by Bloom et al. (1998) on 
August 7 UT. In the case of the GRB 990123 host galaxy 
our BVRcIc observations are the latest. We consider 
that the continuum was expected with minimum change 
due to the fading of the OT. 



^ This paper gives more accurate BVRcIc values with the UT dates 
(corresponding to Table |I|) of observations for the host galaxy of GRB 
980703 unlike the preliminary BVRcIc values given in GCN notice 
#147. 



3. Cosmological models 

The estimate of intrinsic physical parameters of extra- 
galactic objects with redshifts approaching 1 depends 
considerably on the adopted cosmological model. The 
standard Friedmann model contains three parameters: 
Hubble constant Hq, matter density parameter 0,^, and 
cosmological constant parameter Recent studies of 
the Hubble constant put it within the range 50 - 70 km 
s^^ Mpc^^ (see Theureau et al. 1997). In this paper we 
adopt Hq = 60 km s^^ Mpc^^ . The values of Qm and 
Q,\ are observationally less constrained. The recent work 
on the m — z test with supernovae of type la at redshifts 
up to 1 by Garnavich et al. (1998) make it imperative 
to consider in addition to the standard inflationary model 
also an empty universe with a cosmological constant. For 
a review of modem cosmological models and the neces- 
sary mathematical relations, see Baryshev et al. (1994). 
Here we use three Friedmann models which conveniently 
limit reasonable possibilities: 

Ho = 60 km s~i Mpc^^ flm = IMa = (A) 
i?o = 60 km s-i Mpc-\ Jl^ = 0, J^a = (B) 

i^o = 60 km s-i Mpc-\ n,^ = 0, J7a = 1 (C) 

For these models the relation Qm ± + f^fc = 1 is 
valid, where Qrn = Po8ttG/3H^, = Ac^/SH^, and 

= —kc^/RQHQ. Here p, A, k, and R are density, cos- 
mological constant, curvature constant, and radius of cur- 
vature, respectively, and "0" denotes the present epoch. 

The luminosity distance Rium^ the angular size dis- 
tance Rang and the proper metric distance Rp are con- 
nected by the relation: 

Rlu^ = Rangil + zf = Rp{l + Z) (1) 
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where the proper distances for the adopted models are 
given by 



Rn 



jy 2(z-v^T+i+l) 



z{l+0.5z) 
1+z 



Rh 
^Rh z 



for model A, 
for model B, 
for model C. 



(2) 



Here Rh = c/Hq is the present value of the Hubble 
radius. 

The absolute magnitude M(j) of the source observed 
in filter (i) can be calculated from the magnitude-redshift 
relation 

M(,) = m(i) - K^,){z) - 5log{Rium/Mpc) - 25 (3) 

where m^j) is the observed magnitude of the object in 
the photometric band system (i) and K(^i){z) is the K- 
correction at redshift z, calculated from the rest-frame 
spectral energy distribution. 

The linear size I of an object having an angular size 
9 is given by 



I = ORang = ORp/{l + z) 



(4) 



The K-correction in Eq. (3) can be calculated from the 
standard formula (Oke & Sandage 1968): 



= 2.51og(l + z) + 



+ 2.5 loe 



/o°^F,%(A)dA 



(5) 



In this formula F\ is the rest-frame spectral energy dis- 
tribution, is the sensitivity function for the filter (i). 
For our photometric system we used sensitivity functions 
for BVRJc filters from Bessel M. (1990). 

4. Estimation of the K-corrections 

The first spectral observations of the OT of GRB 980703 
were obtained with the Keck-II 10-m telescope on UT 
1998 July 07.6 and 19.6 (Djorgovski et al. 1998). Several 
strong emission lines were detected. These were [Oil], 
US, H7, up and [OIII] with redshift z = 0.9662 ± 
0.0002. In addition, in the blue part of the spectrum some 
absorbtion features (Fell and Mgll, Mgl - absorbtion sys- 
tems) with z = 0.9656 ± 0.0006 were detected. In the 
case of the OT of GRB 990123 absorbtion lines were de- 
tected only (Kelson et al. 1999, Hjorth et al. 1999) with 
Zabs = 1.6004. The HST image reveals that the opti- 
cal transient is offset by 0".67 from an extended object 
(Bloom et al. 1999). This galaxy is most likely to be a 



host galaxy of GRB 990123 and source of the absorbtion 
lines of metals at redshift z = 1.6004. 

In Figure [I] we compared the BVRcIc broad band 
spectra of the host galaxies of GRB 980703 and 
GRB 990123 with the spectral energy distributions of 
normal galaxies of different Hubble types (Pence, 1976). 
Table ^presents fluxes of the host galaxies. Here we used 
the zero-points from Fukugita et al. (1995) for our photo- 
metric bands. For the GRB 980703 the host galaxy fluxes 
are presented according to dereddened magnitudes. For 
the GRB 990123 within our magnitude errors the Galac- 
tic reddening is negligible. The central Xobs = Ae// 
for our photometric system are equal correspondingly 
to: Ab = 4448A, Xy = 5505A, \r = 6588A and 
Xj = 8O6OA, FWHM are equal to: AA^ = 1008A, 
AAv = 827A, A A/?, = 1568A, A A/ = 1542A, respec- 
tively. 

Using equation ^ and the spectral energy distribu- 
tion of the Im Hubble-type galaxies we estimated the 
value of the K-correction for the magnitudes of the host 
galaxies of GRB 980703 and GRB 990123 according to 
z = 0.966 and z = 1.6, respectively. The estimated val- 
ues of the K-correction in the B-hnnA are Kb = 0.68 
for the GRB 980703 host galaxy and Kb = 0.88 for 
the GRB 990123 host galaxy. In this case, the abso- 
lute magnitudes from equation ^ for the host galaxy of 
GRB 980703 are: Mb,,,, = -21.29 for model (A), 
Msr^st = -21.81 for model (B) and Mb,,,, = -22.42 
for model (C). For the host galaxy of GRB 990123 in the 
same way absolute magnitudes are: Mb,,,, = —20.95 
for model (A), Mb,,,, = -21.77 for model (B) and 
Mb,,,, = -22.57 for model (C). 

However, we consider it to be not quite correct to 
compare our broad band spectra to normal Hubble types 
of galaxies. Obviously the starburst activity may drasti- 
caly change the spectral distribution towards the ultravi- 
olet part of the spectrum. According to this considera- 
tion, we compared our BVRcIc spectra to the averaged 
spectral energy distribution of starburst galaxies from 
Connoly et al. (1995). Figure ^ and Figure ^ demostrate 
a comparision of the starburst averaged spectral energy 
distributions to the BVRcIc broad band spectra of the 
host galaxies GRB 980703 and GRB 990123, respec- 
tively. The spectra of sturburst were grouped according 
to increasing values of the color excess E{B — V): from 
SI, with E{B -V) = 0.05 to S6, with E{B -V) = 0.7 
(Connoly et al. 1995). Using relation for (Balmer op- 
tical depth) from Calzetti et al. (1994) we derived the 
values of color excess for individual starburst galaxies. It 
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Figure 1: A comparision of the GRB 980703 and GRB 990123 host galaxy broad band rest-frame (z = 0.966 and 
z = 1.6004 respectively) spectra log Fx = log Fx^obs + C to averaged continuum spectra of galaxies of different 
Hubble types. The spectra were shifted by some arbitrary constants for the best fits. FWHM of each filter for X^jf 
with the account for z are denoted by dashed horizontal lines with bars. 



Table 2: Fluxes of the host galaxies 



Host 


Band 


log FluXx,obs 
erg 

.9-fm2 A 


log FluXy^obs 
erg 
s-crri-^Hz 


GRB 980703 


B 


-17.468 ±0.048 


-28.656 ± 0.048 




V 


-17.508 ±0.040 


-28.509 ± 0.040 




Rc 


-17.588 ±0.032 


-28.440 ± 0.032 




Ic 


-17.823 ±0.072 


-28.483 ± 0.072 


GRB 990123 


B 


-18.168 ±0.064 


-29.356 ± 0.064 




V 


-18.231 ±0.052 


-29.233 ±0.052 




Rc 


-18.456 ±0.056 


-29.308 ± 0.056 




Ic 


-18.579 ±0.120 


-29.247 ±0.120 



are E{B -V) < 0.10 for SI, 0.11 < E{B -V) < 0.21 
for S2, 0.25 < E{B - V) < 0.35 for S3, 0.39 < 
E{B -V)< 0.50 for S4, 0.51 < E{B -V) < 0.60 for 
S5 and 0.61 < E{B -V) < 0.70 for S6 (see Table 3 in 
Calzetti et al. 1994). In Figure |2| and Figure ^ the spectra 
of the SI and S2 type galaxies was averaged with a lOA 
window. 

In this case, the calculations of the K-corrections 
from equation m yield: Kb = —0.01 for the host galaxy 
of GRB 980703 and Kb = 0.13 for the host galaxy 
of GRB 990123. Then, the absolute magnitudes for the 
host galaxy of GRB980703 are: Mb^..^ = -20.60 
for model (A), Ms^^st = -21.12 for model (B) and 
Mb,,,^ = -21.73 for model (C). For the host of 
GRB990123 Mb are: Mb,,,^ = -20.20 for model (A), 



Msrest = -21.02 for model (B) and Mb_, = -21.82 
for model (C). 

5. Estimations of star-forming rate 

We have roughly estimated also SFR using the con- 
tinuum luminosity at Xj-est = 2800A (see Madau et 
al. 1998). In the calculations we assumed the cosmolog- 
ical models described above. 

For the host galaxy of GRB 980703 the effective 
wavelength of the V band for z = 0.966 corresponds 
to 2800A in the rest frame. Using the flux in the V band 
we estimated SFR for the host galaxy of GRB980703: 

SFRsalpeter ~ 15 ± 2 Mq yr-\ SFRscalo ~ 

23 ± 2 Mq yr-1 for model (A); SFRgaipeter ~ 
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Figure 2: A comparision of the GRB 980703 host galaxy broad band rest-frame (z = 0.966j spectrum to spectrum of 
S2-galaxies log Fx = log F\^s2 + C (see Connoly et al., 1995). FWHM of each filter for Ae// with the account for z 
are denoted by dashed horizontal lines with bars. 
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Figure 3: A comparision of the GRB 990123 host galaxy broad band rest-frame (z = 1.6) spectrum to spectrum of 
Sl-galaxies log F\ = log -Fa,5i + C (see Connoly et al., 1995). FWHM of each filter for X^ff with the account for z 
are denoted by dashed horizontal lines with bars. 



24 ± 2 Mo 2/r-\ SFRscaio ~ 37 ± 4 Mq yv'^ 
for model (B); SFRsaipeter ~ 43 ± 4 yr-\ 
SFRscalo ~ 66 lb 6 Mq yr~^ for model (C), where 
the index of Salpeter and Scalo denotes the Salpeter and 
Scalo initial mass function (IMF) (Madau et al. 1998). 

To estimate SFR of the host galaxy of GRB 990123, 
we used the interpolated value of the flux at the wave- 
length 2800A in the rest frame between Rc and Ic band. 
We assumed 

logF ! = -29.28±0.13 . 

^ i/,280oA s-cm?Hz 

The calculation yields: SFRsaipeter ~ 8 ± 2 Mq yr~^. 



SFRscalo ~ 12 ± 4 Mq yr-i for model (A); 

SFRsaipeter ^ 17 ± 6 Mq yr-\ SFRscalo « 25 ± 

9 Mq yr-^ for model (B); SFRgaipeter ~ 34 ± 
12 Mq yr-\ SFRscalo « 54 ± 17 Mq yr'^ for 
model (C). 

Of course, these stimates are the lower limit to SFR 
because our calculations were performed without any 
galaxy rest-frame extinction correction. Moreover, un- 
certaines of our results are estimated formally from the 
errors of fluxes. 
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6. Discussion 

Observations were carried out a significant time after the 
gamma-ray bursts. For the host galaxy of GRB 980703 it 
was about 20 days after the gamma-ray burst, and for the 
GRB 990123 host galaxy about half a year. This allows 
us to consider that the contribution of the optical tran- 
sient is very small and we observed light only of the host 
galaxies. 

As a discussion, it should be noted that there are 
uncertainties in the estimates of the absolute magni- 
tudes of the host galaxies due to the K-correction. In 
the case of the host galaxy of GRB 990123 this un- 
certainty is about l™, and in the case of the host of 
GRB 980703 it is about 0™.7. However, we consider 
that a more correct result is the value of Mb according 
to the comparision to the starburst averaged spectra — 
Msrest = -20.60, -21.12, -21.73 for the host galaxy 
of GRB 980703 in (A), (B), (C) cosmological models re- 
spectively, and MBrest = -20.20,-21.02,-21.82 for 
the host galaxy of GRB 990123. Moreover, our BVRJc 
broad band spectra are better fitted by the S 1 and S2 spec- 
tral energy distribution. To compare our results to the re- 
sults of Bloom et al. (1998,1999) we assume a cosmol- 
ogy model with Hq = 65 km Mpc~^, Qm = 0.2 
and = 0. For the GRB 980703 host galaxy we used 
the value of the luminosity distance from Bloom et al. 
(1998), dL = 1.92 X 10^8 cm. In the case of GRB 990123 
host galaxy we used the value of the luminosity distance 
from Bloom et al. (1999), cIl = 3.7 x 10^^ cm. Cal- 
culations from equation |3| yield: Ms^^^t = —20.8 and 
Mb,,,, = -20.62 for the GRB 980703 and GRB 990123 
host galaxies respectively, while the values of Bloom 
et al. are -20.2 and -20.0 for the GRB 980703 and 
GRB 990123 host galaxies, respectively. Note that the 
estimates of absolute magnitudes of Bloom et al. (1998, 
1999) for the both host galaxies are performed in another 
way without the K-correction by fitting to the spectrum 
in the case of the GRB 980703 host galaxy and by inter- 
polating between the observed the STIS and the K-band 
data points using a power law in the case of GRB 990123 
host galaxy. 

It should be noted that our estimates of the star- 
forming rate are higher than the estimates of Djorgovski 
et al. (1998) and Bloom et al. (1999). It is interesting that 
the for GRB 980703 host galaxy our flux at A = 2800A is 
~ 3.1/iJy and is matching the value of the flux on July 7 
from Djorgovski et al. (1998). Probably, the OT contribu- 
tion was already negligible in the V band on July 7. In the 



case of GRB 990123 host galaxy our value of the flux at 
A = 2800A was estimated by interpolation between ob- 
served points. However, the values of Bloom et al. (1999) 
are 0.17/iJy (for /3 = 0) and 0.21/xJy (for /3 = -0.8) 
while our ones are 0.52~^Q\^fiJy. This discrepancy can 
be explain as follow. Our estimate of the flux was per- 
formed by interpolation between the Rc (Xef / = 6588A) 
and Ic (Xef / = 8O6OA) bands, while the values of Bloom 
et al. were interpolated with power law between STIS 
(approximately V band, Ae// = 5505A) point and K 
band (Ae// = 2.195;um). Moreover, Bloom et al. mea- 
sured the flux of the host galaxy by masking sets of pixel 
dominated by the OT (Bloom et al. 1999) because obser- 
vations in the K band were carried out on epoch 9 and 10 
Febrary UT, 17-18 days after gamma-ray burst, (Bloom 
et al. 1999) when contribution of the OT was not negli- 
gible, while our observations was performed about a half 
year after gamma-ray burst occured. Obviously, the our 
estimate is more exact than that of Bloom et al. 
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